Selmer's multiplicative algorithm 




Abstract 



The behavior of the multiplicative acceleration of Selmer's algo- 
rithm is widely unknown and no general result on convergence has 
been detected yet. Solely for its 2-dimensional, periodic expansions 
exist some results on convergence and approximation due to Fritz 
Schweiger [Sch05]. In this paper we show that periodic expansions of 
any dimension do in fact converge and that the coordinates of the limit 
points are rational functions of the largest eigenvalue of the periodicity 
matrix. 

1 Introduction 

There are at least two ways to approach a theory of multidimensional con- 
tinued fractions (MCFs). One is more of a geometric nature, while another 
concentrates on multidimensional continued fractions which can be described 
by a set of (n + 1) x {n + l)-matrices. The letter set includes amongst others 
the Jacobi- Perron algorithm, as well as algorithms of Brun and Selmer. Each 



to higher dimensions in order to achieve an equivalent of Lagrange's theorem, 
but for cubic or higher roots. However it still remains a challenge, whether 
one of these algorithms applied to an arbitrary ra-tuple x = (ai,...,a n ) 
whereas 1, cvi, . . . , a n are linearly independent over Q and belong to a real 
number field of degree n + 1 eventually becomes periodic. 
Despite this lack of knowledge a lot of discoveries have been made concern- 
ing periodicity and approximation characteristics of multidimensional con- 
tinued fractions. For further information the books of Brentjes [ Bre81] and 




associated with continued fractions 



Schweiger |SchOO] . both named "multidimensional continued fractions" are 
highly recommended. 

1.1 Fibred systems 

In this subsection we first introduce the notion of fibred systems and illustrate 
its characteristics in the context of continued fractions. Thus we are able to 
define multidimensional continued fractions by a set of matrices on such a 
fibred system. 

Definition 1 (fibred system) 

Let B be a set and T: B — ► B be a map. The pair (B, T) is called a fibred 
system if the following conditions are satisfied: 

1. There is a Gnite or countable set I (called the digit set). 

2. There is a map k: B — ► I. Then the sets 

B(i) = fc -1 {i} = {x e B : k(x) = %} 
form a partition of B, hence {J ieI B(i) = B. 

3. The restriction ofT to any B(i) is an injective map. 
(The partition {B(i) : i G 1} is called the time-l-partition) 

Example We consider the continued fraction algorithm. As we know to 
each irrational 16 [0, 1] there corresponds a continued fraction expansion 

x = [a , ai, a 2 ,...] = a -\ , 

a x H ■ 

a 2 + ■■■ 

where we restrict a±, a 2 , . . . to be positive integers. Now we set B := [0, 1[ 
and obtain for i^O the map 

T: B -> B 

x i — ► a(x), a(x) := [x^ 1 \ . 
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For x = however we arrive at a stop and thus get Tx = 0. Further by 
setting I := N we obtain the map i: B —> I defined by 



Then the set / := No is countable, while the sets B(i) form a partition of B. 
Now we set y := Tx = \ — and restrict T to one of the sets B(i),i G I. 
Then for xi, x 2 G B(i) with Tx\ = Tx 2 we immediately obtain x\ = x 2 , since 
a(xi) = a(x 2 ) = i. Hence the pair (B,T) indicates a fibred system. 

With the knowledge of fibred systems we are now able to define multidi- 
mensional continued fractions on the aforementioned systems simply by a 
set of matrices. 

Definition 2 (multidimensional continued fraction) 

The fibred system (B,T) is called (multidimensional) continued fraction if 

1. B is a subset ofW 1 . 

2. For every digit k G / there is an invertible matrix a = a(k) = ((Aij)), 
< i, j < n, such that y = Tx, x G B(k), is given as 



In particular, we are interested in the inverse matrix of such a multidimen- 
sional continued fraction, as hereby we obtain an expression of x via its image 
under the map T. 

Definition 3 

If (B,T) is a (multidimensional) continued fraction, then we denote the in- 
verse matrix of ct(k) by f3(k) = ((By)), < i,j < n. Then we define for 





Aqo + Sj=i AojXj 



s > 1: 



(3(h, . . . , k a ) := (3(h) . . . (3{k s ) = < i, j < n, 



whereas B>- ; = B^. Then y = T s x is equivalent to 



1.2 The concept of cylinders 

In the previous subsection we saw that the digit set / of a fibred system 
causes a partition of the set B into the subsets B(i),i G /, that is 

\jB(x) = B. 

iel 

Simultaneously we didn't demand that T restricted to B(i),i G / is a surjec- 
tive map. Therefore a partition of TB(i),i G / could be of interest as well 
and we eventually arrived at the concept of cylinders. 

Definition 4 (cylinder) 

The cylinder of rank s, defined by the digits i%, . . . , i s is the set 
B(h, . . . , i s ) := B^HT^B^n. . .r\T~^B(i s ) = {x:i 1 (x) = i 1 ,..., i s (x) = i s }. 
Such a cylinder B(i\, . . . ,i s ) is called proper (or full) if T s B{ii, . . . , i s ) = B. 
Proposition 1 

All cylinders of arbitrary rank s, s G N are full if all cylinders of rank 1 are 
full. 

Beweis: Assume that all cylinders of rank 1 are full, then we get 

TB(i) = B 

for all i G I. The rest follows from the Definition Q of cylinders, since for 
(ii, . . . , i s ) G I s we obtain 

TB(ix, ...,i s ) = B(i 2 , . . . ,i s ). 

□ 
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2 Selmer's algorithm 



In connection with Brun's algorithm, Selmer published 1961 a variation of 
its subtractive version called the subtractive algorithm of Selmer (SSA). But 
instead of subtracting the second biggest initial value from the biggest like 
Brun did, he chose to subtract the smallest initial value from the biggest. 
This may at first seem to be a marginal deviation, but it actually implicates 
a fundamental change. As at some point in the expansion one inevitably 
ends up in the absorbing set D := B(n — 1) U5(n). 



2.1 Subtractive version 

Let 

A" +1 := {b = (&oA, . . . A) :b > h >•••>&„> 0}, 
then we define 

ab := (b - 6 n ,&i, . . .AO- 
There is an index i = i(b), < i < n, such that 

nab := (&i, & 2 , . . . , h, b -&„,... , &„) G A n+1 . 

Further let 

B n := {(xi,...,a; n ) : 1 > Xl > ■ ■ ■ > x n > 0}. 
Then with the help of the projection p : A™ +1 — > £> n defined by 



&i 



,&0 &0. 

we obtain the map T : B n — > B n which makes the diagramm 

A n+1 

P 

5 n > B r 



commutative. 
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2.2 Fibred system and absorbing set 



In this subsection we prove that the SSA intrinsically represents a fibred sys- 
tem. Afterwards we shortly point to the absorbing set of Selmer's algorithm. 
As a matter of fact, the restriction of the SSA to this absorbing set coincides 
with an algorithm of Monkemeyer. Hence for further references we point to 
|Mon54| . |Pan08| and jSchOOj . 

Proposition 2 

The Partition 

B(j) := {x E B n : iip^x) = 3} , < j < n, 
makes (B n , T) a fibred system. 

Proof: We set y := Tx and calculate the map T resp. the inverse branches 
< j < n, as follows: 



3 = 



Vk = : Xk , 1 < k < n, x k = Vk , 1 < k < n. 



1 < j < Tl - 1 



Vk = , 1<«<J-1, x\ 



xi yj + y r , 

Vi = — , Xk = Vh ~ l , 2 <k <j, 

Vk = —, 3 + 1 < k < n, x k = — — — , j + 1 < k < n. 
xi Vj + y n 



3 = n 



Xk+l t , . 1 



^1 2/n-l + y 



n 



J/n = , Xfe = ; , 2 < k < n. 

Xl y n -l + Vn 

Therefore we obtain: 
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1. The digit set I := {0, . . . , n} is a countable set. 

2. The map z: B n — > / defined by z: x i— > causes a partition of 
_B n , since evidently = 

3. If we restrict T to any B(j), then we obtain for y = y and < j < n: 



j = As y n = y n it follows that x' n = x" n and thus x = x" . 



1 < j ' < n — 1 As y n = y n and y - = y we obtain x x = a; 1 and hence x = x 



J = n 



As y n = y n and y n _ 1 = y n _ x we get x 1 = x 1 and therefore x = x . 
Consequently the restriction of T to any B(j) is an injective map. 

□ 

Theorem 1 (absorbing set) 

Let D := {x G -B n : x n _i + x n > 1}. Tiien D is an absorbing set, i.e. 

1. TD = D 

2. For almost every x G B n there is an N = N(x), such that T N x G D. 

Proof: A proof of this theorem can be found in |Sch00t p.55] and verifies 
that 

D = B{n- l)n%). 

□ 

2.3 A case of periodicity 

As the question of periodicity is by far the most interesting one, we now give 
an example of a periodic SSA. 

Definition 5 (periodic continued fraction) 

The multidimensional continued fraction of x is called periodic if there are 
numbers m > 0, p > 1 such that T m+P x = T m x. 
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Example (SSA) We consider the tuple x := (x 1 , x 2 ) = (v 7 ? -1,^2-1) 
and apply Selmer's algorithm. Then we obtain 



Tx 



T 2 x 



T 30 x 



T 31 x 



^4-1 $2-1 
2 - $2' 2 - $2 

'3-2^2 $2-1 
#4-1 ' #4-1 



/54- 29^2- 11^4 24^2 + 3^4-35 
^30^4 + 13^2 - 64' 30^4 + 13^2-64 

'27^4-11^2 - 29 24^2 + 3^4-35 
54 - 29^2 - 11^4' 54 - 29^2 - 11#4 



and an easy calculation shows that T 31 x = Tx. Hence the SSA for x := 
(xi,x 2 ) = (#4 — 1, $2 — l) becomes periodic with a preperiod of length 1 
and a period with length 30. 

3 Multiplicative algorithms 

As for Brun's algorithm the multiplicative version really causes an accel- 
eration of expansions, the multiplicative version of Selmer's algorithm in 
general doesn't. In addition none of its cylinders are full. Thus it's more ad- 
equate to denote it as a mere division algorithm. Under these circumstances 
it seems unlikely that the algorithm provides convergent expansions or ap- 
proximations that are competitive in the field of multidimensional continued 
fractions. However, it does. 

3.1 Selmer's division algorithm 

Let A n+1 := {b = (b , b 1 ,...,b n ): b > h > ■ ■ ■ > b n > 0}. Then we define 

5b := (b - kb n ,bi, . . . ,&„), k : = ^ 

Since b n > b — kb n we get n : A™ +1 — > A n+1 defined by 

rr5b := (bi, ... ,b n ,b - kb n ). 
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Now let B n := {(x u ...,x n ): 1 > x x > ■ ■ ■ > x n > 0}. With the help of the 
projection p : A n+1 — > B n defined by 



p(6o,6i, ...,b 
we finally get the bottom map 

S(xi, . . . , x n ) = 
which makes the diagramm 



bo ' " ' ' b 



X2 Xfi 1 fcxn 

5 • • • 5 ) 

X\ X\ X\ 



V 



V 



B n > B n 

S 

commutative. 

Hence the multiplicative version of Selmer's algorithm (MSA) is given by 
S: B n -> B n 

X2 Xfi 1 fcXji 



S(Xi, . . . , Xn) 



J " ■ " J 5 

#1 £1 Xi 



3.2 The fibred system 

In this subsection we simply prove that the MSA represents a fibred system. 

Proposition 3 

The partition 

B(k) : = jx G B n : ^ < x n < jj- J , fc = 1, 2, . . . 
makes (B n , S) a fibred system. 

Proof: If x n = 0, we simply restrict to B" 1 ^ 1 : = {(xi, . . . , x n _i) : 1 > Xi > 
• • • > x n _i > 0}. Then from k := [x" 1 ] we immediately get k G N and thus: 
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1. The digit set / := N is a countable set. 



2. The map i: B n — > / amounts with i: x i— > /c := [a^ 1 ] to a partition of 
5", since evidently U fee/J B(£;) = B n . 

3. If we restrict T to then for y' = y" and thus y n ' = y n ", y n -i = 
y n -i" we obtain x\ = x\ and hence x' = x". Therefore the restriction 
of T to B(k) is an injective map for all fceN. 

□ 

3.3 Cylinders and time-l-partition 

Since k = [x^ 1 ] the pair (B n , S) is a fibred system with cells 

B(k) := ^xeB n : ^ < x n < jj- j , fc = 1, 2, ... . 

These B(k) denote cylinders of rank 1 and indicate in our case convex sets 
with vertices 

,]-),( 



k J \ k + 1 

1 1\ / 1 1 

1,..., T ,T , I,---, 



k' k r\ fc + i' jfe + i 



i i 

k , '" , kj' yk + v"' k + i 

only depending on fc G N. By the way we obtain ^(^j, • • • , ^j) = (1, . . . , 1) 
for all cylinders B(k), fceN. 

In order to achieve more clearness we restrict our attention to the 2-dimensional 
case. Thus we now consider the set 

B 2 = {1 > xi > x 2 > 0}. 

and since 1 > X\ > x 2 und \> x 2 > we obtain the convex set B{k) with 
vertices 



A;/ ' V k + l) 1 \k' k) 1 \k + V k + 1 
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i i i 1 

Figure 1: The time-l-partition of the set B 2 by the 2-dimensional MSA, 
where fceN indicates the associated cylinder B(k). 

These cells B(k) form a partition of the set B 2 that can be easily illustrated 
by means of figure ([I]). 

However none of the cylinders B(k) are full, as they are mapped under S 
onto the convex set with vertices 

(p°)'(^T'*tt)' (1,0MU) ' 

Hence SB(k) C SB(k + 1) for all k G N and additionally SB(k) is not a 
union of cylinders of rank 1. 
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3.4 The matrices 



Notice that if we set y := Sx, then 

1 

Xi = 



2 -/ 



ky„-i + y„ 
Vi-i 



2<i<n. 



ky n -i +y n ' 

Thus, according to the definition of a multidimensional continued fraction, 
the associated (n + 1) x (n+ 1)— matrices of the MSA are given by 

/0 ... k l\ 
1 



and in dimension n = 3 by 



1 





v 



/o 





k 


l\ 


1 














1 













1 





We see immediately, that expanding the determinant along the last column 
leads to 

detp(k) = (-l) n detl = ±1. (1) 
Now we define the matrices P^(k±, . . . , k s ) in common notation as 



B 



^ s \h, ...,k s ):— (3{ki) . . . (3(k s ) 



(s-n+l) 

o 



n(s-i) D 

n n 



(s) 



B 



(s-n) 



R (s-n+l) 
i J->n 



Hence for s > we obtain the relation 

B \ s+l) = k s+l Bt n+1) + B\ s ~ n \ i = 0, . . . , n, 

whereas (3^ denotes the unit matrix. If we set y = S s x and ki = k(S t ~ 1 ), 
1 < i < s, then we find that 



Xi 



(s-n+2) 



j (s-n) 



, i = l,...,n. (2) 
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3.5 Periodic expansions for Selmer's division algorithm 



In this subsection we eventually prove the convergence of the periodic MSA 
and quote a simple example of periodicity. Since this proof allows us to 
apply a variety of other theorems to Selmer's multiplicative algorithm, we 
subsequently mention a few of them. However for completeness and further 
information we refer to [SchOO]. 

3.5.1 Weak convergence of the periodic MSA 
Definition 6 (weak convergence) 

The multidimensional continued fraction is weakly convergent if for every 

x e B 

( b[ s J b$ 

s^oo \ d( s ) ' ' d( s ) 
\ n 00 ^00 

If for arbitrary g with < g < n we get 




s^oo 



B (sy 

R W ' " ' ' R («) 



X, 



for every x G B we call the multidimensional continued fraction uniformly 
weakly convergent. 

Theorem 2 

Assume that the algorithm of x = (x\, . . . ,x n ) eventually becomes periodic 
with period length p. Then 

= X. 

Thus the periodic, multiplicative algorithm of Selmer is weakly convergent 
and even uniformly weakly convergent. 

Proof: Clearly we can assume that the expansion is purely periodic with 

period length p. 

Let 

M:=(3^(h,...,k p ) 




/ r>(p-n+l) td(p) r>(p-n)\ 

I -Dq • • • -Dn -Dn -Dn \ 



\Bir~ n+1) ... B { r 1} b^ B { r n) j 
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Note that 



M 



\-t>n ■ ■ ■ -Dn -£>n -£>n 



Then the characteristic polynomial of M is given as 

Xiviit) := det (tt — M) = t n+1 - b n t n bit - b , 

whereas bo = (— l) n_1 det (M). Furthermore we denote its eigenvalues by 
p , Pi, ■ ■ ■ , p n - Assume that the matrix M has entries equal to 0, then with 
the help of Lemma ([I]) we get m(n), such that M m is a positive matrix for 
all m > m(n) and we could continue with the period length mp. Thus we 
can assume, that M is a positive matrix and the period length remains p. 
Now the Perron- Frobenius theorem |Sch00t Seite 112] allows us to arrange 
the eigenvalues in the way 

Po > \pi\ >•••> \Pn\, Po > 1, 

whereas po is a simple and positive root of Xm(^)- Further with the help of 
the famous Caley-Hamilton Theorem we obtain 

M n+l - b n M n biM -b = 

and by multiplying this with M k f3(ki, . . . , kj) we see that for the entries of 
the matrices the relations 

B ((k+n+l)p+j) _ ^ ^ B ((k+n)p+j) biB {i - k+1)P+j) - b B ikP+j) = 

hold for < i < n and < j < p. Now we apply Theorem 41 of [SchOOj. 
Seite 114] and therefore obtain the general solution 



B? P+3) = d(i,j) P k Q 

+ hoihM + bn(i,j) (i) ft 1 + ■■■ + bi mi (hj) [ mi k _ 1 ) P f mi+1 
+ ... 

^k j_ u fo ^\ f \ „k-l , it /• -\ ( k \ n k-m a +l 



+ b s0 (i,j)p k s + b sl (i,j){ 1 ) p k s +... + b ams (i,j) [ m ^_ 1 jPs 

(3) 
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whereas m 1; . . . ,m s are the multiplicities of the roots pi,...,p s and thus 
1 + mi + • • • + m s — n + 1. If the start values B^ kp+ ^ are given for all k = 
0, . . . , n, then the solution sequence (B^ kp+:jS> ), k > 1 is uniquely determined. 
Now we will consider the terms in equation ^ more precisely. Clearly from 
Po > \pi\ > • • • > \Pn\ an d Po > 1 follows 

Po > |pM > 1 < 7 < s, < fi < m 1 - 1. 

Furthermore for all 7, 1 < 7 < s and e > the relations 

'k + l\ fk\ k+1 t . . 

1 < p < m 7 — 1 



p / \P/ & + 1 — p ' 
(l + £ )*+i = (l +£ )fc(l + £ ) 

hold for all G N and quite evidently follows 

k -\- 1 k^oo 



k + 1- fx 
(l + e)> 1. 



1, 1 < P < W2 7 — 1 



Hence there exists a fc', such that for all k > k' (e) the inequality < 
(1 +e) holds and consequently with increasing k, k > k (s) the term (1 

grows faster than y^j • Thus there exists a fc(e) > fc'(e), such that for all 

k > k(e) the inequality 

^ < (1 + e:) fe V // with 1 < // < m 7 - 1. 

is valid. Since p > 1, there is a e > 0, such that p = 1 + e. Consequently 
there is also a k*(e) > k(e), such that for all k > k*(e) the inequality 



Po>(„)|pM ; 0<p<m,-l. 



holds. So we derive from equation (|3| for d(i,j) 7^ the limit 



lim 1 = d(i,j), (4) 
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as n, n G N is a finite integer and b lfl represents a constant for all 7, p with 
1<7<s,0</!< m 7 , respectively. 
Now we consider the recursion relations 



B 



(s+n) 



B 



(s+rt+l) 



Wib! s+1) +b! !) 



It follows, that if #f s) > pg then b\ s +j) > pg for all j G N, j > n • (n - 1). 
Since d(i,j) doesn't depend on /c and n is finite, it follows that if d(i,j) 7^ 
for some j, then d(i,j) 7^ for all j, < j < p. 
We know that for the trace of M k the relation 



B, 



(fcp-n+1) 




(5) 



holds and since this equation ^ is valid for arbitrarily big k and n is finite, 
there is at least one summand B^ p ~ n+l \ . . . , B^fi, Bn , for whose related 



d(i,j) follows d(i,j) 7^ by equation (J3J). Hence for this i we obtain d(i,j) 7^ 
for all j, < j < p. 

Now we consider the relation for Xi in equation ([2]). As y = T s x G B{k s+ i) := 
\x G -B n : r 1 , : < x n < t~ — \ we get 



1 



< 



B 



(s-n+l) 



< A; s+ i + 1, 



k s+ i + 1 " 5^ n+1) + ■ ■ ■ + B { S) + Sj s - n) 
and by defining k* := max(fci, . . . , k p ) + 1 we eventually obtain 

1 



n 



R (s-n+l) R (s) R (s-n) 

< g ! , B > / 7i ' s < fc* 



Bt n+1) + ■■■ + B$> + B£ 



(s-n) 



n. 



(6) 



Since equation ^ is bounded as well for arbitrarily big s, but d(i,j) 7^ is 
valid for < j < p and a fixed i G {0, . . . , n}, we obtain d(i,j) 7^ for all 
< z < n, < j < p. 
The equation 

/1\ /l\ 



M 



x% 



\X n J 



A 
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shows that for a positive eigenvalue A of the matrix M the relation 

B ^-n + l) + XlB (kp-n+2) + + Xn _ lB (kp) + XnB fP--) = X k . (7) 

holds. But since equation ([7]) is valid for arbitrarily big k and d(0, j) 7^ for 
all < j < p, we get A = p . Thus by equation ^ we derive for T p x = x 
the relation 

Xl = p -* ( B f p ~n+i) + XlB \ kp - n+2) + ■■■ + x n ^B\ kp) + x n Bf p ~ n) ) l<i<n. 



(8) 

As this equation (j8j) holds for arbitrarily big k and d(i,j) exists for all < 
% < n, < j < p, the limit 

Xl = Km p - fc + XljB f p ~ n+2) + ■■■ + x n ^Bl kp) + x n Bf p ~ n) ) (9) 

exists for all i, 1 < i < n. With the help of equation (|4]) we obtain for all 
i, 1 < i < n the relation 

Xi = d(i, p — n + 1) + Xid(i,p — n + 2) + ■ ■ ■ + x^xd^i, 0) + x n d(i,p — n). (10) 

Now we consider the equation j3 2kp+ i = M k f3 kp+ i and see that for the entries 
of the matrices the relations 

B 2kp+j _ B (kp-n+l) j^(kp+j) _|_ _ _ _ _|_ B (kp) g(kp+j) _|_ j^(kp-n) B (kp+j) 

hold for 1 < i < n and < j < p—1. Further we derive Pq > \p k \ > ■ • ■ > 
from p > \pi\ > ■ ■ ■ > \p n \ and hence 

d(i,j) = d(i,p — n + l)d(0, \-d(i, 0)d(n — 1, j) + d(i,p — n)d(n, j). (11) 

follows. 

Now we set 

»> = §§)' (12) 

in equation (10) and obtain equation (11) as a result. By equation (jlj we 
know, that the determinant of M k is unequal to zero. Thus equation ( [To] ) 
represents a system of n equations in n variables and Xi is uniquely deter- 
mined by equation (12). Since this is valid for all j, < j < p and d(i,j) 
doesn't depend on k we get 

□ 
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Lemma 1 

Let M be a (n + 1) x (n + 1)— matrix defined by 



M := 



/O ... k l\ 
1 ... 



1 0) 



whereas k e N represents a positive integer. Then there is a p(n) G N, sucii 
that M p is a positive matrix for all p > p(n). 



Proof: Notice that M — E' + K, whereas 



E' : = 



/0 ... l\ 
1 ... 



K : = 



/0 ... k 0\ 
... 



0/ 



Furthermore (i?') n+1 = 1 represents the unit matrix and K 2 the zero matrix. 
Then apparently the relation 



K(E') 



An-l 



V° ••• o/ 



holds. 

Now we set p = n 2 + 1 and consider M p . Then we get M p = K* + A, whereas 
K* is given by 

(k n ~ l ■■■ k 1 k n \ 



K* 



k n 
1 



n-2 



k 1 

k 

n Un—1 



1 k n k n ~ x ) 
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and A > 0. The shape of K* results from 

n—1 n 

K* = (E') n2+1 + J2 ^2{E'Y[K(E') n ' l ] 3 K(E') n2 ~ jn ' 

j=0 i=0 
n 

,2 , 



i=0 
n 

+ ^( J E , ) l K(E , ) n ~ 1 K( J B , ) ri2 



i=0 

+ ••• 

n 



+ Vf^V KiE'y- 1 . . . K(E') n ~ l K(E') n -\ 
f j \ ✓ 



Hence M p has only posivite entries for all p > n 2 + 1 and therefore (3 P is a 
positive matrix, since k > 1. □ 

3.5.2 Example of a periodic MSA 

Unfortunately the MSA of x :— (v^4, v^2) doesn't become periodic within 
the first 40 steps of expansion. Although we are unaware, whether or not 
periodicity eventually occurs in this expansion, there certainly are periodic 
expansions, some even of period length 1. 

Example (MSA) We consider the tuple x := (x 1 ,x 2 ) = (^p 1 , and 
apply Selmer's multiplicative algorithm. Note that x-i — x\. Then we obtain 



Tx 



>/5-l' y/h-l 
and by multiplying each fraction with ^|+* we obtain 

^"^ 2 ' 2 

Hence the MSA for x := (x 1 ,x 2 ) = (^y^, ^i^) becomes periodic with a 
period of length 1. 



19 



3.5.3 Some general results on the periodic MCF 

Since we proved in Theorem ^ that the MSA is uniformly weakly convergent 
we can show that the coordinates of the limit points are rational functions of 
the largest eigenvalue of the periodicity matrix. Secondary we can apply a 
theorem on approximation properties of MCFs which can be traced back to 
Perron |Per07j . Their general forms for multidimensional continued fractions 
can be found in |Sch001 p. 115-119, 157-162]. 



eventually becomes periodic 



Theorem 3 

Assume that the algorithm of 
with period length p. Then rational functions in p , whereas 

Po denotes the largest eigenvalue of the characteristic polynomial of the peri- 
odicity matrix (3^. Therefore xi,...,x n belong to a number field of degree 
< n + 1. 

Proof: Clearly we can assume that the expansion is purely periodic with 
period length p. Due to Theorem ^ we know, that 

/1\ /l\ 



M 



Xi 



Xi 



and thus get 



(kp—n+l) 




X\By 



(fcp-n+2) r>(kp) | r> 

"T ' " " ~r X n —±±JQ + X n D y 



)(kp-n) 
'0 



Po- 



Therefore we can calculate X\,...,x n as rational functions in po from the 
equation 



M 



\%n J 



PO 



□ 



Below we imply by writing a ^ b with a,6 e K that there is a constant 
C G N, such that a < Cb is valid and accordingly by writing a ^> b with 
a,i G K that there is a constant C G N, such that the inequality a > Cb 
holds. 
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Theorem 4 

For n-dimensional, periodic algorithms we obtain for arbitrary e > the 
estimates 

|^ P9) x 4 - J B i (ra) |«| Pl (l + e)^ i = l,...,n. 
But there is at least one pair such that the inequality 

holds for infinitely many values of g. 

Proof: Due to Theorem ^ we get for all i, 1 < i < n the relation 
B^Xi - B^ P9) = k(i, 0)pg 

+ m 10 (i, 0)p? + m u {i, 0) (A pf 1 + ■■■ + m lmi (i, 0) ( 9 _ ^ p? 

+ ... 



+ m s0 (t, 0)pf + m sl (i, 0) I J 1 pf" 1 + • • • + m sms (z, 0) I y _ J p 



and with the help of equation (12) it follows, that k(i,0) = for all i, 1 < 
i < n. Thus the relation reduces to 



^ - = rn 10 (i, 0)p? + m u (i, 0) I \\ pf 1 + • • • + m lmi (i, 0) [ *_ J p? 



'/•;/! ... 77'/"/' mi() /,' r»\„f . .... m i 'M „.'/-! . /; n\ / y 1 ,,.</-»< i + i 

+ . . . 



+m s0 (z, 0)pf + m al (i, 0) (J J pf- 1 + ■ • • + m sms (z, 0) 1 J p 

By Theorem ^ we know, that for any e > there exists a such that 
for all g > g(e) the inequality 



<(l + e) 9 , l<p<m 7 -l, 1<7<s. 



holds. Since for the eigenvalues |pi| > • • • > \p s \ holds, Ipi] 9 > ■ ■ • > |p s | 9 is 
valid as well, where s < n. Hence there is a constant c(x,e), such that for 
all g > g(s,c) the inequality 

\BMxt-BW\Kc\n0. + e)\°, i = l,...,n. 
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holds. This proves the first part of our theorem. If all roots pi, . . . , p& with 
Px — ■ ■ • — Pd are simple, then 



\B x { - B i | 



n. 



The second part of the proof can be found in [SchOOl Theorem 59]. 



□ 
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